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Abstract 

The non-diagonal correlators of vector and scalar currents are considered at three- 
loop order in QCD. The full mass dependence is computed in the case where one 
of the quarks is massless and the other one carries the mass M. As applications 
we consider the single-top-quark production via the process qq — > tb and the decay 
rate of a charged Higgs into hadrons. In both cases the computed NLO corrections 
are shown to be numerically much less important than the leading ones. 
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1 Introduction 



In high energy experiments where the center-of-mass energy, s/s, is much larger than the 
mass of the quarks, M, the latter can often be neglected or an expansion in Mjsfs is 
sufficient to describe the experimental data. However, there are situations where the full 
dependence on M and yfs is required. One can, e.g., think on the high precision which 
meanwhile has been reached at LEP (CERN), SLC (SLAC) or TEVATRON (Fermilab) 
or on situations where the center-of-mass energy is of the same order of magnitude as the 
quark masses. In particular for threshold phenomena the masses are important and the 
full dependence is desirable. 

A variety of important observables can be described by the correlators of two currents 
with different tensorial structure. If the coupling of the currents to quarks is diagonal 
quantities like e + e~ annihilation into hadrons and the decay of the Z boson are covered 
by the vector and axial-vector current correlators. Total decay rates of neutral CP even 
or CP odd Higgs bosons can be obtained with the help of the scalar and pseudo-scalar 
current densities, respectively. For these cases the full mass dependence at order a 2 has 
been computed in JTJ for the non-singlet and in J| for the singlet correlators. 

In this work we deal with three-loop correlators of currents which couple to two dif- 
ferent quarks flavours: a massless one and one carrying the mass M. 

The non-diagonal vector and axial-vector correlators describe properties connected to 
the W boson. The absorptive part of these correlators is directly related to the decay 
width of a (highly virtual) W boson into quark pairs and gluons. Of particular interest 
in this connection is the single-top-quark production via the process qq — > tb H |J. 
This is because the process directly probes the matrix element \Vtb\ - one of the few 
parameters of the Standard model not yet experimentally measured. The imaginary part 
of the transversal W boson polarization function constitutes a gauge invariant and finite 
contribution of 0{a 2 ) to the process. 

As an application of the (pseudo-)scalar current correlator we consider the decay of a 
charged Higgs boson which occurs in extensions of the Standard Model. The corrections 
provided in this paper describe the total hadronic decay rate of such a boson (coupled to 
a massive and a massless quark). 

2 Definitions 

In the vector case the polarization function is defined through 



(-?V + ^) Tl v (q 2 ) + q,q u ni(q 2 ) = i J dx e^{0\Tj;(x)fJ (0)\0), (1) 



with j v = V'lT/iV^- Only the transversal part U v (q 2 ) will be considered in the following. 
The definition of the scalar polarization function reads 




(2) 
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with j s = (m(fi) I 'M)%pi%p2 where m(/x) is the MS quark mass and M is its pole mass. 
Throughout this paper we consider anti-commuting 75. This is justified because only 
non-singlet diagrams contribute (-01 is supposed to be different from ^2). Thus, within 
perturbation theory, the axial-vector and pseudo-scalar correlators coincide with the vec- 
tor and scalar ones, respectively. 

We consider only the case where one of the quarks is massive. Thus we will identify 
tpi with q, the massless quark, and ip 2 with Q which is supposed to be a heavy quark of 
mass M. Furthermore it is convenient to introduce the dimensionless variable 

z = h> (3) 

where M refers to the pole mass. For the overall normalization of H 5 (q 2 ) (5 = v, s) we 
adopt the QED-like conditions IT 5 (0) = 0. 

The physical observables R(s) are related to IT(g 2 ) through 

nV = a + ie)l, (4) 



R v {s) = 12vr Im 
R s (s) = 8vr Im 
where the use of the variables 



nV = s + *e) , (5) 



M 1-x 2 
x = v = — — - , 6 

y/S 1 + X 2 

turns out to be useful to describe the high energy and threshold region, respectively. 
The expansion of II s in terms of ot s reads (5 = v, s) 

a ^V) ^ u (i),s + I «f__V ) N 
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u s = U (0),s + *_JW ^(i)^ + ^Atl + 0(a 3 s ) . (7) 



It is convenient to decompose the three-loop term according to the colour structure 

IT^ = C 2 F uf/ + C A C^P/ + C F T ni uf£ S + C F Tuf^ s , (8) 

("2) S 

where analogous formulae hold for R(s). In Eq. (H\ U F p corresponds to the abelian part 
already present in QED whereas the non-abelian structure is contained in 11^ . The 
remaining two structures correspond to the fermionic contributions where n\ counts the 
number of massless quarks and nj = n\ + 1 is the total number of active quark flavours. 



3 Calculation 

A complete analytical computation of Il <5 (g 2 ) at three-loop order or its imaginary part is 
currently not feasible. The method we use for the computation of the diagrams is based on 
conformal mapping and Pade approximation [||, |!| [I], H . It allows for the computation of 
a semi-numerical approximation for Il <s (g 2 ). The aim is the reconstruction of the function 
II 5 (g 2 ) from the knowledge of some moments for z — > and z — > —00 and additional 
partial information about the behaviour of R(s) close to the threshold, i.e. for s — ► M 2 . 
Thus let us in the following briefly discuss the different kinematical regions. 
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Expansions for small and large external momentum 

The number of diagrams which contribute at one-, two- and three- loop order is relatively 
small and amounts to one, three and 31, respectively. Nevertheless we used GEFICOM 0] for 
the automatic computation. GEFICOM uses QGRAF || for the generation of the diagrams. 
In case an asymptotic expansion has to be applied LMP |{J or EXP jlO| are used for the 



generation of the sub-diagrams. The occurring vacuum diagrams are passed to MATAD pi 
and the massless propagator type diagrams are evaluated with MINCER [[[ 
on the automatic computation of Feynman diagrams can be found in |~3] 



More details 



We have been able to compute seven terms for small and eight terms for large external 
momentum q both for the vector and scalar correlators of Eqs. (]!]) and (|2|), respectively. 
This means expansion coefficients up to order z 6 , respectively, 1/z 7 are available for the 
Pade procedure. The analytical results are rather bulky and will be published elsewhere. 



Threshold behaviour 

The Born results for the vector and scalar spectral functions are proportional to v 2 in the 
limit v — > as can be seen for the exact results 



R (0), 
R (0), 
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(9) 
(10) 



Also the corrections of order a s to R" and R s are known [lj, [15], [L£], [T7|, [L8| and a 
similar threshold behaviour is observed: the spectral functions vanish like v 2 (modulo 
powers of lnt> ). This is valid in every order in a s as follows from Heavy Quark Effective 
Theory (HQET) ||19|| . Actually the latter can be used to obtain the leading threshold 
behaviour of R v (s) and R s (s) at 0(a 2 ) from the corresponding correlators in HQET. In 
particular, the renormalization group equation in the effective theory is used to get the 
leading logarithmic behaviour at order a 2 s . Afterwards the decoupling relation between 
the currents in the full and effective theory fl20| , [H| , 22| is exploited to get the information 
about R v (s) and R s (s). For jx 2 = M 2 our results read 



R 



v,thr 



N c v 



+ 



»i6 + a -°"'Wg,(^ + 126-9h,2-9h» 

7T V 2 



a 



(n/) (M) 



7T 



cl 



'FF 



147 27 \ 27 
30C 2 + — In 2 In v + — In 2 v 



+ C A C F [c FA + 



423 33, \, 33, 2 

19( 2 H In 2 }lnv H In 2 v 

8 S 2 J 4 



+ C F Tm {c v FL + 
133 



39 



+ 8C 2 -61n2 lnw-31n 2 w 



+ C F T 



- 10c 



IF 



3 



,(D,v 




,(l),s 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

v 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

v 



Figure 1: R^' v (s) and R^ ,s (s) function of v. 
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(12) 



with = tt 2 /6. The constants Cp F ,c^ A and c FL (5 = s,v) are unknown. Note, however, 
that the leading term for the colour structure CpT is completely determined. In order 
to incorporate the available threshold information into out method one has to perform 
an analytical continuation of the expressions in Eqs. (|TTD and fll2p. Taking the logarith- 
mic parts of Eqs. (|Tl~D and ( |T2"D one obtains the quadratic and cubic logarithms for the 
polarization functions. They can be incorporated into the Pade procedure n. 



4 Results and Applications 

In a first step we want to discuss the results for the spectral functions R" and R s which 
essentially corresponds to the direct outcome of the Pade method. For the results we 
present in the following, those Pade approximants are chosen which contain for their 
construction at least terms of order z 5 and 1/z 5 in the small and large momentum region, 
respectively. Furthermore we demand that the difference of the degree in the numerator 
and denominator is less or equal to two. In the numerical evaluations we specify fi 2 = M 2 . 
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Figure 2: R F p(s), RfJ"(s), R FI f(s) and R F h(s) as a function of v. 

The analytical formulae which result from the semi-numerical Pade proce- 
dure are quite long. Thus we refrain from listing them explicitly. Instead, 
a typical representative for each colour structure can be found under the URL 
http : / /www-ttp . physik . uni-karlsruhe . de/Progdata/ ttp00-25. 

At order a s the exact results are known. Thus it is instructive to make a comparison 
with the outcome of our semi-numerical method. In Fig. [I] roughly 15 Pade approximants 
are plotted together with the exact curves. Note that even for small values of v, where only 
the logarithmic parts of Eqs. ([TT|) and ([12]) are incorporated in our method, no difference 
both between the individual Pade expression and the exact result is visible. This is quite 
promising for the 0{a 2 s ) results. The corresponding curves are shown in Fig. |2] and Fig. |3| 
for the vector and scalar case, respectively. Again of the order of 15 Pade approximants 
are plotted and no difference is visible. We should note that due to our choice /x 2 = M 2 
in the high energy limit all curves tend to ±oo except Rpp which reaches a constant for 
v — ► 1. In particular, Rp\ s and Rp^ s tend to +oo, respectively, — oo. This is not visible 
in the figures as the corresponding range in v where the turn over takes place is beyond 
the resolution in the plots. 

In order to get even more confidence in the results we plot in Fig. |] Rp F v (s) and 
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Figure 3: R F p(s), Rfa(s), R f £ s (s) and R FI f(s) as a function of v. 

R^Pp S (s) as a function of x. Expressed in this variable the high energy region gets more 
spread. Thus a sensible comparison with the expansion terms can be performed. The 
dashed curves correspond to the analytical expressions obtained via asymptotic expansion 
containing the terms up to order 1/z 7 . Excellent agreement with the semi- numerical 
results is observed up to x ~ 0.5 for the vector case and x ~ 0.7 for the scalar case which 
corresponds to v ~ 0.60 and v ~ 0.34, respectively. This is a strong consistency check for 
the Pade method. 

Single-top-quark production 

As a first application we want to discuss the single-top-quark production via the process 
qq — > tb. Some sample diagrams are plotted in Fig. [|. 

The corrections of order a s to the (total) single-top-quark production rate are quite 
large. They amount to about 54% and 50% for Tevatron and LHC energies, respec- 
tively P3J, where 18%, respectively, 17% arise from the final state corrections. This calls 
for a complete 0(a 2 s ) calculation. In this letter we want to do a first step and consider 
the leading term in the large- N c expansion. 
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Figure 5: Sample diagrams contributing to the process qq — > tb. The wavy and loopy 
lines represent W bosons and gluons, respectively. 

At 0(a s ) there is no interference between the initial and final state radiation. Thus 
it is possible to write the differential cross section is factorized form 

da ( - +1,y\ < x*r* , v ^^w(q 2 ,M 2 ,M 2 ) 

— (pp^tb + X) = a( PP ^W +X) <q2 _ m2 , (13) 

where Tlw corresponds to the transversal part of the W boson correlator and is connected 
to the vector correlator of Eq. ([]]) through 

IW) = V2G F M^\V tb \ 2 q 2 U v (q 2 ). (14) 



Due to interference diagrams between the initial and final state (c.f. Fig. |^(f)) Eq. (|13D 
does not hold at a 2 . However, those contributions are suppressed by at least a factor 
1/N 2 in large N c limit as compared to the diagrams in Fig. |](c)-(e). The latter, together 
with the contributions in Fig. |5|(a) and (b), are covered by Eqs. ([13]) and (14). 
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At order a 2 there are also diagrams like the one in Fig. |5](g) which appear for the 
first time. In principle they lead to the same final state as the diagram in Fig. [|(e). 
However, one has to note that the W boson generating the top and bottom quark is 
radiated from a light quark flavour. Furthermore, compared to the diagram in Fig. |5](c) 
these contributions are only suppressed by a factor 1/N C and not by 1/N 2 like the diagram 
in Fig. |(f). 

Thus, if we restrict ourselves to the leading term in 1/N C it is possible to use the results 
for R v obtained above in combination with Eq. (|l^) to perform a theoretical analysis of 
the order a 2 to the single-top-quark production in the large- N c limit. In order to obtain 
the total cross section the corresponding parton distribution functions would be needed 
to the same order. 

The production cross-section of the virtual W* boson is identical to that of the Drell- 
Yan process qq — > ev e . The latter is known to 0(a 2 ) from Ref. [24]. Thus we can take the 



proper ratios to make predictions in the large- iV c limit at NNLO free from any dependence 
on parton distribution functions. As an example, we consider 

cicr 



^(pp^W*^tb) _ Im[H 46 (g 2 )] 



(pp — > W* — > ez/ e ) ImTUg 2 )] 

= N c \V tb \ 2 R v (s). (15) 

To get an impression of the numerical significance we plot in Fig. the LO, NLO and 
NNLO result of R v (s) in the range yfs = 200 . . .400 GeV. For the numerical values we 
choose M t = 175 GeV and a s (M z ) = 0.118. Whereas the 0(a s ) corrections are significant 
there is only a moderate contribution from the order a 2 terms. In the range in q 2 shown 
in Fig. |6| they are below 1% of the Born result. Note that the NNLO correction to the 
Drell-Yan process are also small and amount to at most a few percent (see e.g. |25||). 
Thus, in case there is no kinematical magnification for the diagrams in Fig. [5](f) and (g) 
we can conclude that the radiative corrections to the single-top-quark production via the 
process qq — > tb are well under control. 



Decay rate of a charged Higgs 

Theories beyond the SM are usually characterized by an enlarged Higgs sector and may 
allow for different quantum numbers of the Higgs bosons. For example, one of the most 
appealing extensions of the SM, the Minimal Supersymmetric Standard Model (MSSM), 
contains two complex iso-doublets with opposite hyper-charge (see e.g. resulting in 

five mass eigenstates of (pseudo-)scalar physical Higgs fields: two neutral CP-even (H° 
and h°), one neutral CP-odd (A) and two charged (H^) Higgs bosons. 

Let us consider a generic charged Higgs boson coupled to fermions through 

Ch+dv = (V2G F ) 1/2 H + J H+ , (16) 
where the corresponding quark current is given by 

J H+ = !^U[a(l- l5 ) + b(l+ l5 )}D. (17) 
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Figure 6: LO (dotted), NLO (dashed) and NNLO (solid) results of R v (s). 

Here U and D represent generic up- and down-type quarks, respectively, with MS masses 
mjj and thd = 0. We only consider H + as for H~ the formulae are analogues. The 
parameters a and b are model dependent and are left unspecified. 

The decay rate of the Higgs boson H + into quarks and gluons can be written in the 
form 

T(H + ^UD) = V2G F M H+ lm \li H {M 2 H+ )\ , (18) 

where Mu is the pole quark mass and Uh (q 2 ) is given by 

q 2 U H (q 2 ) = Jdxe^ x (TJ H + (x)J H -(0)) = (a 2 + b 2 )q 2 U s (q 2 ) , (19) 

Thus, we arrive at the following expression for the hadronic decay rate of the charged 
Higgs boson 

Y(H + ^UD) = ^LM H+ {a 2 + b 2 )R s {M 2 H+ ). (20) 



37T 

In Fig. R s (Mj i+ ) is plotted at LO, NLO and NNLO. Again it turns out that the radiative 
corrections are well under control as order a 2 terms contribute at most of the order of 
1%. 
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Figure 7: LO (dotted), NLO (dashed) and NNLO (solid) results of R s (s), M — M t — 
175 GeV. 

5 Conclusions 

We have computed the non-diagonal correlators of vector and scalar currents at order a 2 
in the case of the currents composed from a massless and a massive quark field. Our main 
theoretical tools have been the asymptotic expansions in the large and small momentum 
region as well as some information about the threshold behaviour provided by HQET. 

The obtained results for the case of the (axial-)vector correlator constitute an impor- 
tant ingredient of the full calculation of next-to-next-to-leading 0(a 2 ) correction to the 
single-top-quark production via the process qq — > tb. Such a calculation is required for an 
accurate extraction of the matrix element \Vtb\ from experiment. 

We have used the results for the (pseudo-) scalar correlator in order to obtain the 0(a 2 ) 
decay rate of a charged (pseudo-) scalar Higgs boson to hadrons. 
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